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An improved leading order fragmentation functions set of kaon is proposed based on the experi- 
mental data. We compare it with currently available sets, and use it to calculate high-pr Kj-K ratios 
in relativistic proton-proton collisions. A prediction of Kj-K ratios in relativistic nucleus- nucleus col- 
lisions (d -I- Au and Au -\- Au) at RHIC i/s = 200 GeV is given. 
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An enhancement of the abundance of strange quarks 
(s and s) is expected if a quark-gluon plasma (QGP) is 
formed in a relativistic nuclear collision, compared to a 
collision without QGP formation^. This asymmetry in 
the flavor composition of QGP may be reflected in the 
final particle composition. This is therefore of interest to 
systematically study the production of strange mesons, 
in particular kaons, in relativistic proton-proton (pp) and 
proton-nucleus ipA) collisions. Furthermore, it has been 
suggested 3] that the measurable Kj-K ratios are sensitive 
to the initial density of the QGP in relativistic nucleus- 
nucleus i^AE) collisions. Henceforth, in order to provide 
a basis for the treatment of kaon production in nuclear 
collisions at the energies of the Relativistic Heavy Ion 
Collider (RHIC), it is important and necessary to first 
develop a good description of kaon production in col- 
lisions and relying on that, we can then give a reliable 
studying of hard (high-pT^) Kj-K probe for nuclear colli- 
sions. 

In the framework of parton model of perturbative QCD 
(pQCD), the invariant cross section for inclusive kaon 
production in a pp collision gets involved with fragmenta- 
tion functions (FFs), D}^/a.{z, Q^), through the factoriza- 
tion theorem The value of Dh/aiZjQ"^) corresponds 
to the probability for the parton a produced at short 
distance ~ l/Q to form a jet that includes the hadron 
h carrying the fraction z of the longitudinal momentum 
of parton a. It is a nonperturbative object and hence 
can not be derived from QCD. During the last few years, 
the experimental progress at the CERN and the SLAG 
has motivated people to generate leading-order (LO) and 
next-to-leading-order (NLO) set ofn^, K^, and p/p FFs 
set ^IMI^i which can be utilized in pQCD calculations. 
However, the available FFs sets for kaon are still not sat- 
isfied. For example, in Fig^ we show the comparison 
between Kretzer's FFs (K) and that of Kniehl, Kramer, 
and Potter (KKP), for the fragmentation process of u 
quark to + K~. We find that the differences between 
KKP-FFs and K-FFs are significant, especially at large-z 
region. 

On the other hand, the experimental data of kaon 
production in pp collisions in the energy range 20 GeV 
< < 40 GeV0, shows that K+/7r+ 1/2 at 
high-pT region, which reminds us that /D^^ 1/2. 
While interesting enough, it has been implied by Feyn- 
man et al. that, when z ^ 1 (z is the longitudinal mo- 



mentum fraction of a parton), D^* / ^ (1 — P)/f3, 
in which /3 is a constant parameter in Feynman-Field 
parameterization^ 0|. In current letter, based on the 
KKP-FFs of pion and experimental data, we propose the 
following improved LO FFs (Z) ZJf of kaon, 



a, KKP ) 
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in which D'^ kkp ^^e corresponding LO KKP-FFs of 
pion and parton a can be the u{u), d{d), s(s) quark or 
gluon (g). 

To deal with charged kaons, i.e., K'^ and K~ , we now 
consider the scheme to separate the above FFs for kaon. 
A natural way to do is as follows: we assume, say, for u 
valence quark of K'^ (us) , 
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so that we have 



D'^ — — D'" 

^u,Z — 2^^KKP J 



(2) 



(3) 



which means the contribution of u valence quark to 
K~ (us) is vanished. 

Starting with eqn-s.|(21) and we compare our im- 
proved Z-FFs with KKP-FFs and K-FFs in FigUl We 
can see that the Z-FFs agrees with K-FFs in the most 
of moderate z region, but disagrees with it at large-z re- 
gion (with z > 0.7), where in stead tends to agree with 
KKP-FFs. 

Based on this Z-FFs set, the LO invariant cross sec- 
tion of kaon production in pp collisions is described in 
the pQCD-improved parton model on the basis of the 
factorization theorem as a convolution llOl: 
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J dxa dxb dkxa dkxb dz^ 

abed 



X9a/pO^Ta) fa/p{Xa,Q^) gb/pO^Tb) fb/p{xb, Q^) 



-sd{s + i + u), (4) 



x—^{ab^cd) — 

dt TTZ^ 



where fa/p{x,Q^) and fb/p{x,Q'^) are the parton distri- 
bution functions (PDFs) for the colliding partons a and 
h in the interacting protons as functions of momentum 
fraction x, at scale Q, da/dt is the hard scattering cross 
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FIG. 1: The comparisons of improved FFs (full line) with 
KKP-FFs (dashed line) and K-FFs (dot-dashed line). In 
top panel we show the fragmentation function of u quark to 
-'rK' , +^ . In the bottom panel, the ratios between 
different FFs sets, Z/KKP (fuU line), K/KKP (dashed line), 
and Z/K (dot-dashed line), are shown, as a function of z at 
fixed scale = 25 GeV^ 



FIG. 2: The comparison of experimental (symbols) and 
pQCD-calculated (curves) K/-k ratios: /-k^ (top panel) 
and K~ (bottom panel), as a function of pt at -^i = 
19.4, 23.7, 27.4, and 38.8 GeV, respectively, in pp collisions. 
In top panel, we only show K'^ /ir'^ ratios at = 27.4 (full 
line) and 38.8 (dashed line) GeV. The kaon production is cal- 
culated based on improved Z-FFs. The experimental data of 
K/tt ratios are from Refs.0, 



section of the partonic subprocess ab — s- cd, and the FFs, 
^K/c(^c,Q'^) gives the probability for parton c to frag- 
ment into kaon with momentum fraction Zc at scale Q' . 
The superscript "Z" is used here for our improved Z- 
FFs. We use the convention that the parton-level Man- 
delstam variables are written with a 'hat' (like t above). 
The scales are fixed in the present work as Q = pt/2 
and Q' — pt/2zc. The function ga/pO^r) is the intrinsic 
transverse momentum distributions of parton a inside nu- 
cleon. We follow the phenomenological approach in the 
present work, taking g(k.T) to be a Gaussian: 

,(k.) . ^e-^/^'^^) . (5) 

Here (fc^) is the 2-dimensionaI width of the fc^ distri- 
bution and it is related to the magnitude of the average 
transverse momentum of one parton as (k^) = 4(A;T)^/7r. 

In this investigation we use LO partonic cross sections, 
together with LO PDFs (GRV) |l3. The Monte-Carlo 
integrals are carried out by the VEGAS-routine . 

We analysed experimental data of inclusive K produc- 



tion in pp collisions in the center-of-mass energy range 
20 GeV < < 40 GeV|l 0. We find a good agree- 
ment between our model and experimental data on both 

/ir^ and /n^ . In our analysis, the intrinsic trans- 
verse momentum distribution of parton, (k'^) = 1.7, 2.0, 
2.2, and 2.2 GeV^, are used for both kaon and pion at 

= 19.4, 23.7, 27.4, and 38.8 GeV, respectively [13 . 

High-pT kaon production data begun to emerge from 
RHIC^at center-of-mass energies ^/s = 130 and 200 
GeV 0, 0|. This provides us with the opportunity 
to predict kaon production in d -I- Au and Au + Au col- 
lisions at higher energies. For nucleus-nucleus collisions, 
the Glauber model[16j is the standard description. The 
LO inclusive kaon production in AB collisions can be 
written in the following formlism: 

^Y.J d^bd\TA{r)TB{\h-r\) j dx, 

abed 

X dXbdkxadkTb dZc ga/AO^Ta)fa/A{Xa, Q^) 
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FIG. 3: The predicted K/n ratios for d + Au and Au + Au 
collisions at RHIC energy ^ = 200 GeV (0 - 5% centrality) 
with HIJING shadowing, but without jet quenching. The 
values of {k'^)pp = 1.7 GeV^, Vm = 4, and C = 0.4 GeV^ are 
fixed 1131 in the calculation. 



X gb/si^Tb) fb/B{xb,Q^) — —{ab -> cd) 



l'2 



dt 



sS{s + t + u). 



(6) 



The predicton on K/tt ratios at RHIC energy, y/s = 
200 GeV, for both d + Au and Au + Au collisions, are 
shown in Fig|3| 

In above calculations, we include parton niultiscatter- 
ing due to the effect of nuclear medium. The extra con- 
tribution to the width can be related to the number of 
nucleon-nucleon (NN) collisions in the medium as fol- 
lows: 



PP 



(7) 



GeV'^ at = 200 GeV [ij. The function hpA{b) de- 
scribes the number of effective NN collisions at impact 
parameter h which impart an average transverse momen- 
tum squared C and it can be written in terms of the 
number of collisions suffered by the incoming proton in 
the target nucleus, i^a(&) = o-NNtA{b), where aNN is 
the inelastic NN cross section. In this work, we take 



the 'saturated' prescription 0,0]., i.e., we assume that 
only one associated NN collision is responsible for the 
enhancement, so h''J^^{h) was equated to a maximum 
value Vrn, whenever VAyb) > I'm = 2. We fix the value of 
= 4 and an associated value of C = 0.4 GeV^ [l3| . 
we also approximately include shadowing effect (the 
PDFs are modified due to nuclear environment) [isl Il9| 
and the isospin asymmetry of heavy nuclei into the nu- 
clear PDFs. We use the average nuclear dependence and 
the scale-independe nt p arameterization of the HIJING 
shadowing function 



nt p 



fa/Aix.Q"^) = Ra/A{x) 



jfa/p{x,Q^) 



+ Ra/Aix) 



fa/n{x,Q') 



(8) 



where fa/ni^i Q^) is the PDF for the neutron. 

In conclusion, an improved LO FFs set of kaon pro- 
duction is proposed in present letter, based on the ex- 
perimental data. Relying on this improved FFs set, we 
calculate the high-pT K/n ratios iupp collisions in center- 
of-mass energy range 20 GeV ^ ^/s ^ 40 GeV, and we 
find our results in pp collisions are in good agreement 
with experiment data. We also predict K/tt ratios for 
d -I- Au and Au + Au collisions at y/s = 200 GeV to pro- 
vide background information for further comparison with 
the high-pT kaon production data of RHIC. The future 
systematic investigation of kaon production with this im- 
proved FFs set in pp, pA and AB collisions will be very 
interesting and be reported in other place. For AB colli- 
sions with the consideration of jet quenching, the current 
FFs set of kaon will give very important insight about the 
existence and properties of the QGP. 

The author thank to Dr. Fai and Dr. Kretzer for useful 
comments and stimulating discussions. 
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